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We estimate the strong coupling constant from the perturbative expansion of the plaquette. The scale is set by 
the 2S-1S and IP-IS splittings in bottomonium which are computed in NRQCD on dynamical gauge configurations 
with rif = 2 degenerate Wilson quarks at intermediate masses. We have increased the statistics of our spectrum 
calculation in order to reliably extrapolate in the sea-quark mass. We find a value of a-g^{mz) = 0.1118(26) 
which is somewhat lower than previous estimates within NRQCD. 



1. Introduction 

In recent years lattice calculations of the strong 
coupling constant have achieved a remarkable 
precision. This is mainly due to the availabil- 
ity of gauge configurations that partially include 
the effects of vacuum polarisation. A promising 
approach has been pioneered by the NRQCD Col- 
laboration [jjj. They examine expectation values 
of small Wilson loops that are very easily and ac- 
curately computed on the lattice and for which 
perturbative expansions are known to 0{a 2 s ). 
The preferred choice to fix the lattice scale, a -1 in 
this context are the lowest radial and orbital split- 
tings in the T system since they do not require 
a precise tuning of the b-quark mass, nib and are 
expected to exhibit little dependence on the light 
quark masses. In this short note we present a de- 
termination of as along the lines of Ref.|j]] that 
is able to clarify two important issues: first, we 
use dynamical Wilson quarks instead of the stag- 
gered formulation, so that one may estimate the 
impact of different discretisation schemes. Sec- 
ond, we calculate the dependence of T levels on 
the dynamical quark mass, m q and thus reduce 
the systematic error from simulating overly heavy 
sea- quarks. 
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configurations 
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0.1560 
0.83 

206 
824 



0.1570 0.1575 

0.76 0.69 

192 203 

2496 2436 



Table 1 

Simulation set-up. 

We rely on the recipe of tadpole improvement us- 
ing uq from the mean-link in Landau gauge to ac- 
count for the bulk of loop contributions and oth- 
erwise stick to tree-level couplings. Propagators 
are evaluated on SESAM's gauge configurations 
with two degenerate flavours of Wilson quarks 
0. We exploit each configuration more than 
once starting at different sources, see Table @. 
We did not tune the bare b-quark mass but kept 
ami, = 1.7 throughout the simulation. This value 
reproduces the correct T mass in the quenched 
approximation and it turns out to be adequate 
in the full theory, too, leading to kinetic masses 
m ki „(T) = 9.97(28), 9.63(24), 9.68(27) GeV for 
k = 0.1560,0.1570,0.1575, respectively. Com- 
pared to the results presented in || , we have sig- 
nificantly increased our statistics except for the 
heaviest sea-quark mass. 



2. Simulation 

We employ NRQCD at next-to-leading or- 
der, i.e. spin-independent operators of 0(m,bV 2 ), 
0(rribV 1 ) and spin-dependent terms of 0(m&u 4 ) 
and 0(mbV 6 ) have been included into the action. 



3. Bottomonium Spectroscopy 

The results quoted here are taken from simul- 
taneous multi-exponential fits to correlators with 
a local sink and a source that is smeared with a 
carefully chosen potential model wave function. 
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splitting aAEo c aAE(m s /3) 

2 L S - 1 L S Q 0.209(21) L2(7) 0.229(10) 

2 3 5_i - 1 3 ^ 0.209(18) 1.1(7) 0.226(9) 

l 3 F-l 3 5i 0.163(9) 0.4(3) 0.170(5) 

2 1 P 1 - ^Pi 0.152(24) 1.8(7) 0.181(15) 

Table 2 

Results of the extrapolation in the sea-quark 
mass. 3 P is the spin-averaged triplct-P state. 

a-^GeV] R^. 

V - T x - T 

7^ = 0,/? = 6.0 2.29(11) 2.68(9) 1.49(10) 

k = 0.1560 2.12(16) 2.38(8) 1.44(10) 

k = 0.1570 2.34(9) 2.50(8) 1.37(7) 

k= 0.1575 2.43(8) 2.57(7) 1.35(6) 

m s /3 2.49(10) 2.59(7) 1.33(6) 

Table 3 

Determination of the lattice spacing from the 
2 3 S'i - l 3 ^ and l 3 P - l 3 Sx splittings. Their 
ratio, Rsp is to he compared to the experimental 
value of 1.28. 

These correlators have the cleanest signals and 
lead to very stable fits. Following Q we extrapo- 
late the energy splittings linearly in the sea-quark 
mass, 

a AE = a AEq + c am q (1) 

to am e ff = am s /3 ~ 0.0159 which is somewhat 
below our lightest quark mass. Fit parameters are 
listed in Table ^ together with the energy split- 
tings at m ff. We find the dependence on the 
dynamical quark mass to be much smaller than 
in the light hadron sector. Results for lattice 
spacings are summarised in Table |3[ We use the 
average spacing from 2 3 Si — l 3 Si and 1 3 P — l 3 Si 
(at m q = m e ff in the unquenched case) to convert 
our results to physical units. Note that the aver- 
age a -1 for (3 = 5.6, nf = 2 at am e g agrees well 
with the quenched one at (3 = 6.0, so that we can 
directly compare both theories. Figure |l| collects 
our spectrum results for two and zero flavours. It 
is obvious that the gross level structure computed 
on quenched configurations disagrees with exper- 
iment whereas the n/ = 2 results are in much 
closer agreement. This is also evident from Table 
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Figure 1. Bottomonium spectrum - radial level 
splittings. The 3 Si ground state is hxed to 
the experimental T mass. Open symbols denote 
quenched results, Riled symbols nt = 2 results 
at m q = m s /3. Solid lines mark experimental 
values, dashed lines the positions of the spin- 
averaged 3 Pj states. 

H which shows that the lattice spacings from the 
2S-1S and the IP-IS splittings do not agree in the 
quenched theory whereas they coincide when two 
dynamical quarks are switched on. 

4. Plaquette Coupling 

We adopt the definition of the plaquette cou- 
pling given in [EI, 

-lnWu = -ln(-3?TrD) = 

O 

4tt /3.41\ r s 
—a P ( — J [1 - (1.1855 + 0.0249 n f )a P }(2) 

which equals the expansion in ay but shifts trun- 
cation errors to the conversion into standard cou- 
plings later on. In Table ^ we summarise the 
couplings cip obtained from Eq. (|J) as well as the 
scales determined from the IP-IS (% — S) and 2S- 
1S (T' — T) splittings. In the unquenched case we 
quote values for both am q = am s /3 and the chiral 
limit to estimate the systematic error connected 
to the finite sea-quark mass. Subsequently these 
couplings are evolved to a common scale using 
the universal two-loop (3 function. The plaquette 
couplings in the quenched and unquenched the- 
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M q 


-lnWn 




Ml [GcV] 








X-T T' - T 


0.0159 


0.5570 


0.1677 


8.84(44) 8.48(46) 





0.5546 


0.1668 


9.16(83) 9.09(101) 


oo 


0.5214 


0.1518 


9.13(50) 7.82(54) 



Table 4 

Results for ap (3.41/a) extracted from the mea- 
sured plaquette values. 

ory can now be extrapolated to the active num- 
ber of light quark flavours which is expected to 
be rif = 3 in the case of the low-lying bb bound 
states. Guided by the perturbative evolution, we 
extrapolate Op 1 linearly in rif, Figure |2| Obvi- 
ously the mismatch between ap-values obtained 
from different splittings in the quenched approx- 
imation disappears, once the dynamical quarks 
are switched on. 

5. Discussion 

To make the connection with the MS'-scheme 
one invokes 



MS 
(n 



(Q) 



o 



5/6 



Q 



i 



-o 



(«/) 



C 2 (n f )c 



where C2 = 0.96 for nt = 0. Following Ref. 
JTJ, we use this value also for n/ = 2 and take 
the whole size of the quenched two-loop contri- 
bution as an error estimate for the conversion. 
The couplings are evolved to the Z mass apply- 
ing third-order perturbative evolution. We per- 
form the matching at the heavy-quark thresholds 
m c — 1.3 GeV and nib = 4.1 GeV. We thus arrive 
at our final results 



"-MS 



{m z ) 



0.1118 (16) (5) (20) x-T 
0.1124 (18)(15)(20) T'-T 



(3) 



The first error includes the statistical uncertainty 
and the systematic errors due to the influence of 
relativistic corrections and changes in uo on the 
lattice scale. The second error quantifies the sea- 
quark mass dependence and the last one the trun- 
cation error in the conversion. 
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Figure 2. The plaquette coupling ap as a func- 
tion of the number of dynamical flavours. The 
triangles result from an extrapolation in the in- 
verse flavour number. 

We cannot confirm the staggered result fij] but 
obtain a value of a-jjg{mz) which is smaller by 
two standard deviations. A closer inspection re- 
veals that the effect can be traced back to the 
difference between staggered and Wilson plaque- 
tte values which is larger than anticipated from 
perturbation theory, while the scales in both sim- 
ulations are comparable. There is some contro- 
versy about the proper way to determine the lat- 
tice scale in an unquenched simulation. Whereas 
SESAM's strategy is to set a -1 in the chiral limit, 
one may also think of separate lattice spacings for 
each quark mass. A priori it is conceivable that 
the value of ap would be larger if it were de- 
termined on each subensemble and then extrap- 
olated in m q . However, we observe a very mild 
dependence of the IP — 15 splitting on the quark 
mass, so that different extrapolation methods do 
not produce a big effect. 
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